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Abstract 

Given an n-vertex planar directed graph with real edge lengths and with no negative 
cycles, we show how to compute single-source shortest path distances in the graph in 
0(n log 2 n/ log log n) time with 0(n) space. This is an improvement of a recent time bound 
of 0(nlog 2 n) by Klein et al. 

1 Introduction 

Computing shortest paths in graphs is one of the most fundamental problems in combinatorial 
optimization with a rich history. Classical shortest path algorithms are the Bellman-Ford algo- 
rithm and Dijkstra's algorithm which both find distances from a given source vertex to all other 
vertices in the graph. The Bellman- Ford algorithm works for general graphs and has running 
time 0(mn) where m resp. n is the number of edges resp. vertices of the graph. Dijkstra's al- 
gorithm runs in 0{m + n log n) time when implemented with Fibonacci heaps but it only works 
for graphs with non-negative edge lengths. 

We are interested in the single-source shortest path (SSSP) problem for planar directed 
graphs. There is an optimal O(n) time algorithm for this problem when all edge lengths are 
non- negative [3]. For planar graphs with arbitrary real edge lengths and with no negative 
cycles, Lipton, Rose, and Tarjan [7] gave an 0(n 3 / 2 ) time algorithm. Henzinger, Klein, Rao, and 
Subramanian [3] obtained a (not strongly) polynomial bound of 0(n 4//3 ). Later, Fakcharoenphol 
and Rao [2] showed how to solve the problem in 0(n log 3 n) time and 0(n log n) space. Recently, 
Klein, Mozes, and Weimann [6] presented a linear space 0(n log 2 n) time algorithm. 

In this paper, we improve the result in [6j by exhibiting a linear space algorithm with 
0(n log 2 nj log log n) running time. 

From the observations in [6], our algorithm can be used to solve bipartite planar perfect 
matching, feasible flow, and feasible circulation in planar graphs in 0(n log 2 nj log log n) time. 

The organization of the paper is as follows. In Section [2 we give some definitions and some 
basic results, most of them related to planar graphs. Our algorithm is very similar to that of 
Klein et al. so in Section [31 we give an overview of some of their ideas. We then show how to 
improve the running time in Section HJ Finally, we make some concluding remarks in Section [5j 



2 Definitions and Basic Results 

In the following, G = (V,E) denotes an n-vertex planar directed graph with real edge lengths 
and with no negative cycles. For vertices u, v G V, we let dc(u, «) elU denote the length 
of a shortest path in G from u to v. We extend this notation to subgraphs of G. We will 

*Department of Computer Science, Brown University, shay@cs.brown.edu, 
http : / /www . cs . brown . edu/people/ shay/ 

^Department of Computer Science, University of Copenhagen, koolooz@diku.dk, 
http: //www. diku. dk/hjemmesider/ ansatte/koolooz/ 



assume that G is triangulated such that there is a path of finite length between each ordered 
pair of vertices of G. The new edges added have sufficiently large lengths so that shortest path 
distances in G will not be affected. 

Given a graph H, we let Vh and Eh denote its vertex set and edge set, respectively. For an 
edge e G Eh, we let 1(e) denote the length of e (we omit H in the definition but this should not 
cause any confusion). Let P = U\, . . . , u m be a path in H. We let \P\ = m. For 1 < i < j < m, 
P\ Ui, Uj\ denotes the subpath Ui, . . . , Uj. A path P 1 is said to intersect P if Vp !~l V P < / 0. If 
P' — u m , . . . , u m ' is another path, we define PP' — u\, . . . , u m —\, u m , ■u m +i, . . . , 

Define a region R (of G) to be the subgraph of G induced by a subset of V. In G, the 
vertices of Vr that are adjacent to vertices in V\ Vr are called boundary vertices (of R) and the 
set of boundary vertices of R is called the boundary of R. Vertices of Vr that are not boundary 
vertices of R are called interior vertices (of R) . 

The cycle separator theorem of Miller [8] states that, given an m-vertex plane graph, there 
is a Jordan curve C intersecting 0(y/m) vertices and no edges such that between m/3 and 2m/3 
vertices are enclosed by C. Furthermore, this Jordan curve can be found in linear time. 

Let r G (0, n) be a parameter. Fakcharoenphol and Rao [2] showed how to recursively apply 
the cycle separator theorem such that in 0(n log n) time, (a plane embedding of) G is divided 
into 0(n/r) regions with some nice properties: 

1. each region contains at most r vertices and 0( v / r) boundary vertices, 

2. no two regions share interior vertices, 

3. each region has a boundary contained in 0(1) faces, defined by simple cycles. 

We refer to such a division as an r-division of G. The bounded faces of a region are its holes. 
To simplify the description of our algorithm, we will refer to all vertices of the cycles containing 
the boundary of a region as boundary vertices of that region. Furthermore, we will assume that 
for each region R in an r-division, R is contained in the bounded region defined by one of the 
cycles C in the boundary of R. Clearly, this can always be achieved by adding a new cycle if 
needed. We refer to C as the external face of R. 

For a graph H, a price function is a function p : Vh — ► R- The reduced cost function induced 
by p is the function w p : Eh — ► R, defined by 

w p (u, v) = p(u) + l(u, v) — p(v). 

We say that p is a feasible price function for H if for all e G Eh, w p (e) > 0. 

It is well known that reduced cost functions preserve shortest paths, meaning that we can 
find shortest paths in H by finding shortest paths in H with edge lengths defined by the reduced 
cost function w p . Furthermore, given (j) and the distance in H w.r.t. w p from a u G Vh to a 
v G Vzj, we can extract the original distance in H from u to v in constant time [6]. 

Observe that if p is feasible, Dijkstra's algorithm can be applied to find shortest path 
distances since then w p (e) > for all e G Eh- An example of a feasible price function is 
u i — > dn(s, u) for any s G Vh- This assumes that dn(s, u) < oo for all u G Vh which can always 
be achieved by, say, triangulating H with edges of sufficiently large length so that shortest paths 
in H will not be affected. 

A matrix M = (Mij) is totally monotone if for every i,i',j,f such that i < i', j < j' and 
Mij < Miji, we also have M^j < M^y . Totally monotone matrices were introduced by Aggarwal 
et al. in [9], who gave an algorithm, nicknamed SMAWK, that, given a totally monotone n x m 
matrix M, finds all column minima of M in just 0(n + m) time. A matrix M = (Mij) is convex 
Monge if for every i,i',j,f such that i < i' , j < j', we have My + Mjy > M^/ + M^j. It is 
immediate that if M is convex Monge then it is totally monotone. Thus SMAWK can be used 
to find the column minima of a convex Monge matrix. [6] used a generalization of SMAWK to 
so called falling staircase matrices, due to Klawe and Kleitman [4], that finds all column minima 
in 0(ma(n) + n) time, where a(n) is the inverse Ackerman function. 



3 The Algorithm of Klein et al. 



In this section, we give an overview of the algorithm of [6] before describing our improved 
algorithm in Section [H 

Let s be a vertex of G. To find SSSP distances in G with source s, the algorithm starts by 
applying the cycle separator theorem to G. This gives a Jordan curve C which separates G into 
two subgraphs, Go and Gi, and there are 0(^/n) boundary vertices on G. 

Let r be any of these boundary vertices. The algorithm consists of five stages: 

Recursive call: SSSP distances in G% with source r are computed recursively for i = 0, 1. 

Intra-part boundary distances: The distances in G% between each pair of boundary vertices 
of Gi are computed using the algorithm of [5] for i = 0,1. This stage takes 0(n log n) time. 

Single-source inter-part boundary distances: A variant of Bellman-Ford is used to com- 
pute SSSP distances in G from r to all boundary vertices on G. The algorithm consists of 
0(y / n) iterations and each iteration runs in 0(y / na(n)) time using the algorithm of Klawe and 
Kleitman [3]. This stage therefore runs in 0{na{n)) time. 

Single-source inter-part distances: Distances in the previous stage are used to modify 
G such that all edge lengths are non-negative without changing the shortest paths. Dijkstra's 
algorithm is then used in the modified graph to obtain SSSP distances in G with source r. Total 
running time for this stage is 0{n log n). 

Rerooting single-source distances: A price function is obtained from the computed dis- 
tances from r in G. This price function is feasible for G and Dijkstra's algorithm is applied to 
obtain SSSP distances in G with source s in 0(n log n) time. 

4 Improved Algorithm 

As can be seen above, the last four stages of the algorithm in [6] run in a total of 0(n log n) 
time. Since there are O(logn) recursion levels, the total running time is 0(n log 2 n). We now 
describe how to improve this time bound. 

The idea is to reduce the number of recursion levels by applying the cycle separator theorem 
of Miller not once but several times at each node of the recursion tree. More precisely, for a 
suitable p, we obtain an ?i/p-division of G in 0(n log n) time. For each region Ri in this n/p- 
division, we pick an arbitrary boundary vertex ri and recursively compute SSSP distances in 
Ri with source r^. This is similar to stage one in the original algorithm except that we recurse 
on 0(p) regions instead of just two. 

We will show how all these recursively computed distances can be used to compute SSSP 
distances in G with source s in 0(n log n + npa(n)) additional time. This bound is no better 
than the 0(n log n) bound of the original algorithm but the speed-up comes from the reduced 
number of recursion levels. Since the size of regions is reduced by a factor of at least p for each 
recursion level, the depth of the recursion tree is only 0(logn/logp). It follows that the total 
running time of our algorithm is 



O 




To minimize this expression, we set nlogn = npa{n). Solving this, we get p = log n/a(ra) which 
gives a running time of 0(n log 2 n/ log log n), as requested. 



What remains is to show how to compute SSSP distances in G with source s in 0(n log n + 
npa(n)) = 0(n log n) time, excluding the time for recursive calls. 

So assume that we are given an n/p-division of G and that for each region R, we are given 
SSSP distances in R with some boundary vertex of R as source. Note that the number of regions 
is 0(p) and each region contains at most n/p vertices and 0(y/n/p) boundary vertices. 

We will assume in the following that no region has holes. Then all its boundary vertices are 
cyclically ordered on its external face. We consider holes in Section 14.41 

The remaining part of the algorithm consists of four stages very similar to those in the 
algorithm of Klein et al. We give an overview of them here and describe them in greater detail 
in the subsections below. Each stage takes 0(n log n) time. 

Intra-region boundary distances: For each region R, distances in R between each pair of 
boundary vertices of R are computed. 

Single-source inter-region boundary distances: Distances in G from an arbitrary bound- 
ary vertex r of an arbitrary region to all boundary vertices of all regions are computed. 

Single-source inter-region distances: Using the distances obtained in the previous stage 
to obtain a modified graph, distances in G from r to all vertices of G are computed using 
Dijkstra's algorithm on the modified graph. 

Rerooting single-source distances: Identical to the final stage of the original algorithm. 

4.1 Intra-region Boundary Distances 

Let R be a region. Since R has no holes, we can apply the multiple-source shortest path 
algorithm of [5] to R since we have a feasible price function from the recursively computed 
distances in R. Total time for this is 0(|Vr| log |Vr|) time which is 0{n log n) over all regions. 

4.2 Single-source Inter-region Boundary Distances 

Let r be some boundary vertex of some region. We need to find distances in G from r to all 
boundary vertices of all regions. To do this, we use a variant of Bellman-Ford similar to that 
in stage three of the original algorithm. 

Let 1Z be the set of 0{p) regions, let B C V be the set of boundary vertices over all regions, 
and let b = \B\ = 0{j)\J njp) = O(^frvp). Note that a vertex in B may belong to several regions. 

Pseudocode of the algorithm is shown in Figure HJ Notice the similarity with the algorithm 
in [6] but also an important difference: in [6], each table entry ej[v] is updated only once. Here, 
it may be updated several times in iteration j since more than one region may have v as a 
boundary vertex. For j > 1, the final value of ej[v] will be 

ej[v]= mm{e j - 1 [w]+ d R (w,v)}, (1) 

where B v is the set of boundary vertices of regions having v as boundary vertex. 
To show the correctness of the algorithm, we need the following two lemmas. 

Lemma 1. Let P be a simple r-to-v shortest path in G where v G B. Then P can be decomposed 
into at most b subpaths P = P1P2P3 . . ., where the endpoints of each subpath Pi are boundary 
vertices and Pi is a shortest path in some region oflZ. 

Proof. P is simple so it can use a boundary vertex at most once. There are b boundary vertices 
in total. A path can only enter and leave a region through boundary vertices of that region. □ 



1. initialize vector ej[v] for j = 0, . . . , b and v G B 

2. ej[v] := oo for all v G B and j = 0, . . . , b 

3. eo[r] :=0 

4. for j = 1, . . . , 6 

5. for each region R G 7£ 

6. let C be the cycle defining the boundary of R 

7. ej[v] := min{ej[i>],min w6 y c {ej_i[u;] + dji(w,v)}} for all i> G Vc 

8. := e 6 [u] for all v G B 

Figure 1: Pseudocode for single-source inter-region boundary distances algorithm. 

Lemma 2. After iteration j of the algorithm in FigureUl ej[v] is the length of a shortest path 
in G from r to v that can be decomposed into at most j subpaths P = P1P2P3 . . . Pj, where the 
endpoints of each subpath Pi are boundary vertices and Pi is a shortest path in a region oflZ. 

Proof. The proof is by induction on j. We omit it since it is similar to that in [6j. □ 

It follows from (JXj) and Lemmas [T] and [2] that after b iterations, D[v] holds the distance in G 
from r to v for all v G B. This shows the correctness of our algorithm. 

Line 7 can be executed in 0(\ Vc\a(\ Vc\)) time using ideas from [6J (since all boundary 
vertices of R are cyclically ordered on its external face) and the fact that dji(w,v) has been 
precomputed in the previous stage for all v,w G Vc- Thus, each iteration of lines 4-7 takes 
0(bot{n)) time, giving a total running time for this stage of 0(b 2 a(n)) = 0(npa(n)). Recalling 
that p = logn/a(n), this bound is O(nlogn), as requested. 

4.3 Single-source Inter-region Distances 

To compute distances in G from boundary vertex r to all vertices of G we consider one region 
at a time. So let R be a region. We need to compute distances in G from r to each vertex of R. 

Let R' be the graph obtained from R by adding a new vertex r' and an edge from r' to 
each boundary vertex of R; the length of this edge is equal to the distance in G from r to the 
boundary vertex. Note that da(r,v) = dR>(r',v) for all v G Vr. Also note that R' has 0(|Vr|) 
vertices and edges and can be computed in 0(|Vr|) time, given the distances computed in the 
previous stage. We need to find distances in R' from r' to each vertex of Vr. 

Let tr be the boundary vertex of R for which distances in R from rR to all vertices of R 
have been recursively computed. We define a price function <f> for R' as follows. Let Br be 
the set of boundary vertices of it! and let D = m&x{dR(rR, b) — dc(r, b)\b G Br}. Then for all 
v G Vr', 



Lemma 3. Function (p defined above is a feasible price function for R' . 

Proof. Let e = (it, v) be an edge of R' . By construction, no edges end in r' so v 7^ r' . If u 7^ r' 
then (j){u) + l(e) — <j){v) = dR(rR,u) + l(u,v) — dR(rR,v) > by the triangle inequality so assume 
that u = r' . Then v G Br so 4>(u) + /(e) — <p(v) = D + da(r, v) — dR(rR, v ) > by definition of 



Price function <p can be computed in time linear in the size of R and Lemma [3] implies that 
Dijkstra's algorithm can be applied to compute distances in R' from r' to all vertices of Vr in 
0(|VrJ log \Vr\) time. Over all regions, this is 0(n log n), as requested. 

We omit the description of the last stage where single-source distances are rerooted to 
source s since it is identical to the last stage of the original algorithm. We have shown that 




D. This shows the lemma. 



□ 



all stages run in O(ralogn) time and it follows that the total running time of our algorithm is 
0(n log 2 n/ log log n). It remains to deal with holes in regions. 

4.4 Dealing with Holes 

In Sections 14. II and 14.21 we needed the assumption that no region has holes. In this section, we 
remove this restriction. As mentioned in Section [21 we may assume w.l.o.g. that each region of 
1Z has at most a constant h number of holes. 

Intra-region boundary distances: Let us first show how to compute intra-region boundary 
distances when regions have holes. The reason why it works in 0(n log n) time in Section [4.11 
is that all boundary vertices of each region are on the external face, allowing us to apply the 
multiple-source shortest path algorithm of [5]. 

Now, consider a region R. If we apply [5] to R we get distances from boundary vertices on 
the external face of R to all boundary vertices of R. This is not enough. We also need distances 
from boundary vertices belonging to the holes of R. 

So consider one of the holes of R. We can transform R in linear time such that this hole 
becomes the external face of R. Having done this transformation, we can apply the algorithm 
of [5] to get distances from boundary vertices of this hole to all boundary vertices of R. If we 
repeat this for all holes, we get distances in R between all pairs of boundary vertices of R in time 
0{\V R \ log \V R \ + h\V R \ log \V R \) = 0{\V R \ log |Vr|) time. Thus, the time bound in Section O 
still holds when regions have holes. 

Single-source inter-region boundary distances: What remains is the problem of com- 
puting single-source inter-region boundary distances when regions have holes. Let C be the 
external face of region R. Let Hr be the directed graph having the boundary vertices of R as 
vertices and having an edge (n, v) of length cIr(u, v) between each pair of vertices u and v. 

Line 7 in Figure [1] relaxes all edges in Hr having both endpoints on C. We need to relax 
all edges of Hr. In the following, when we say that we relax edges of R, we really refer to the 
edges of Hr. 

To relax the edges of R, we consider each pair of cycles (C%, C2), where C\ and C2 are C or 
a hole, and we relax all edges starting in C\ and ending in C2. This will cover all edges we need 
to relax. 

Since the number of choices of (C±,C2) is 0(h 2 ) = 0(1), it suffices to show that in a single 

iteration, the time to relax all edges starting in C\ and ending in C2 is 0((| Vc 1 \ + \Vc 2 \) a fflc 1 1 + 

|Vc7 2 |)), with 0(\ Vr\ log I Vr\) preprocessing time. 

We may assume that C\ 7^ C2, for otherwise we can relax edges as described in Section [4. 2i 
In the following, we define graphs, obtained from R, that are needed in our algorithm. It is 

assumed that these graphs are constructed in a preprocessing step. Later, we bound the time 

to construct them. 

We transform R in such a way that C\ is the external face of R and C2 is a hole of R. We 
may assume that there is a shortest path in R between every ordered pair of vertices, say, by 
adding a pair of oppositely directed edges between each consecutive pair of vertices of Cj in some 
simple walk of Cj, i = 1, 2 (if an edge already exists, a new edge is not added). The lengths of 
the new edges are chosen sufficiently large so that shortest paths in R and their lengths will not 
change. Where appropriate, we will regard R as some fixed planar embedding of that region. 

Before proceeding, let us give some additional definitions. We say that an edge e = (u, v) 
with exactly one endpoint on path P emanates right (left) of P if (a) e is directed away from 
P, and (b) e is to the right (left) of P in the direction of P. If e is directed towards P, then we 
say that e enters P from the right (left) if (v,u) emanates right (left) of P. We extend these 
definitions to paths and say, e.g., that a path Q emanates right of path P if there is an edge of 
Q that emanates right of P. 



Figure 2: Region Rp is obtained from P essentially by cutting open at P the "ring" bounded 
by C\ and C2. 

Now, let n G Vc*! and let T be a shortest path tree in R from n to all vertices of Ci- Let 
P be a simple path in T from ri to some leaf T2 £ Vc 2 - Define E resp. E as the set of edges 
that either emanate left resp. right of P or enter P from the left resp. right. 

Now, take a copy Rp of R and remove P and all edges incident to P in Rp. Add two copies, 
P and P, of P to Pp. Connect path P resp. P to the rest of Rp by attaching the edges of E 
resp. E to the path, see Figure[2l If (u,v) S -Br, where (v,u) £ Ep, we add (u,v) to P and P 
in Pp. We extend this construction to paths from C\ to C2 other than P. 

In order to relax edges from boundary vertices of C\ to boundary vertices of C2 in R, the 
first step is to relax edges in Rp, defined in the following. 

Note that a simple, say counter-clockwise, walk (u\ = r±),U2, ■ ■ ■ , i^icM) = r i) °f Ci 

in R starting and ending in r\ corresponds to a simple path P\ = u'i, . . . , uj Ci | +1 in Rp. In the 
following, we identify Ui with u\ for i = 2, . . . , \C\\. Vertex u\ = r\ in R corresponds to two 
vertices in Pp, namely u'i and u\d\+i- We will identify both of these vertices with r\. 

A simple, say clockwise, walk of C2 in R from r2 to T2 corresponds to a simple path P2 = 
Vi, ... ,v'^ c , +l in Pp. We make a similar identification between vertices of C2 and P2. 

In the following, when we e.g. say that we relax all edges in Rp starting in vertices of C\ 
and ending in vertices of C2, we really refer to edges starting in the corresponding vertices of 
Pi and ending in the corresponding vertices of P^. More precisely, suppose we are in iteration 
j. Then relaxing an edge from a u G Vc 1 \ {n} to a v € Vc 2 \ {^2} in Rp means updating 

ej[v] := mm{ej[v],ej-i[u] + d Rp (u',v')}. 

If u = r±, we relax w.r.t. both u[ and and if v = r%, we relax w.r.t. both v[ and wL^. 

We extend these definitions to graphs with a structure similar to Rp. 

As the following lemma shows, relaxing edges in Rp can be done efficiently by exploiting 
the cyclic order of boundary vertices of Rp as we did above for regions with no holes. 

Lemma 4. Relaxing all edges from Vc\ to Vc 2 in Rp can be done in 0(|Vci| + | Vc- 2 1 ) time in 
any iteration of Bellman-Ford. 

Proof. Let paths P\ and P2 in Rp be defined as above. Consider iteration j. Define a |Pi| x IP2I 
matrix A with elements Aki = ej-i[uk] + dp p {u' k ,v[). Observe that relaxing all edges from Vc\ 
to Vc 2 in Rp is equivalent to finding all column-minima of A (compare this to [B]). 

Now, since P1PP2P is a cycle, it follows easily from results of [6] that for 1 < k < k' < |Pl| 
and 1 < I < I' < IP2I, Am + ^k'V > ^-kV + ^k'l- Hence A is Monge (see [15]), so by [Sj, its 
column-minima can be found in OdVc 1 !) + | Vcr 2 1) time. □ 

Unfortunately, relaxing edges between boundary vertices in Rp will not suffice since shortest 
paths in R that cross P are not represented in Pp. To overcome this obstacle we identify two 



n n 

(a) (b) 

Figure 3: (a): The rightmost root-to- leaf simple path P r and the leftmost root-to- leaf simple 
path Pg in T. (b): In the proof of Lemma El if Q first crosses P r from right to left and then 
crosses Pg from right to left then there is a u-to-v shortest path in R that does not cross Pg. 



particular paths P r and Pg such that for any u G C\ , v G C2 there exists a shortest path in R 
that does not cross both P r and Pg. Then, relaxing all edges between boundary vertices once in 
Rp r and once in Rp e suffices to compute shortest path distances in R. We now formalize this 
idea. 

Let us define what we mean when we say that path Q = q±, 52, #3, . . . crosses path P. Let 
onto be the smallest index such that q ou t does not belong to P. We recursively define irn to 
be smallest index greater than ouij_i such that qi ni belongs to P, and outi to be smallest index 
greater than irij such that qi ni does not belong to P. We say that Q crosses P from the right 
(left) with entry vertex Vi n and exit vertex v ou t if (a) Vi n = qi ni and v ou t = q ou u-i for some i > 
and (b) qi ni -i<lini enters P from the right (left) and (c) q uti-iQout emanates left (right) of P. 

A rightmost (leftmost) path P in T is a path such that no other path Q in T emanates 
right (left) of P. Let P r and Pg be the rightmost and leftmost root-to-leaf simple paths in 
T, respectively; see Figure ^ a) . Let v r G C2 and vg G C2 denote the leaves of P r and Pg, 
respectively. 

Lemma 5. For any u G Vq 1 and any v G Vc 2 , there is a simple shortest path in R from u to v 
which does not cross both P r and Pg. 

Proof. Suppose for the sake of contradiction that there are vertices u G Vc± and v G Vc 2 such 
that any simple shortest path in R from u to v crosses both P r and Pg. Let Q be a simple u-to-v 
shortest path in R such that the sum of the number of times it crosses P r and the number of 
times it crosses Pg is minimal. 

Assume that Q crosses P r first. The case where Q crosses Pg first is symmetric. Let Wi n and 
w ou t be the entry and exit vertices of the first crossing, see Figure E^b). If Q first crosses P r 
from left to right, then observe that it must cross Pg at the same vertices. In fact, it must be 
that all root-to-leaf paths in T coincide until w ou t and that Q crosses all of them. In particular, 
Q crosses the root-to-v path in T, which we denote by T v . Since T v does not cross P r , the 
path Q[u, w ou t]T v [w ou t, v] is a shortest path from u to v in R that does not cross P r . But this 
contradicts our assumption above. 

We conclude that Q first crosses P r from right to left. Consider the path S = Q[u, w 0U t]P r [w ut, v r ]. 
We claim that Q does not cross S. To see this, assume the contrary and let w' denote the exit 
point corresponding to the crossing. Since Q is simple, w' ^ Q[u,w ou t\. So w' G P r [w ou t, v r ], 
but then Q[u, w out ]P r [w ou t, w']Q[w' , v] is a shortest path from u to v in R that crosses P r and 
Pg fewer times than Q. But this contradicts the minimality of Q. 

Since Q first crosses P r from right to left and never crosses S, its first crossing with Pg must 
be right-to-left as well, see Figure [3^b). This implies that Q enters all root-to-leaf paths in T 



1. initialize vector ej [v] for j = 0, . . . , b and v G B 

2. ej [v] := oo for all v G B and j = 0, . . . , b 

3. e o [r]:=0 

4. for j = I, . . . ,b 

5. for each region R £ 1Z 

6. for each pair of cycles, Ci and C2, defining the boundary of R 

7. if Ci = C2, relax edges from C\ to C2 as in [6] 

8. else (assume C\ is external and that dR Pr and g?r p ^ have been precomputed) 

9. ej[v] := min{ej[v],mm we v c { e j-i[ w } + dn Pr (w,v)}} for all v G 

10. e^hv] := min{ej[v],min^ e y c {ej_i[w] + dp p (w,v)}} for all u G Vc 2 

11. D[v] := e b [v] for all v G 5 



Figure 4: Pseudocode for the Bellman-Ford variant that handles regions with holes. 

before (not strictly before) it enters Pi. In particular, Q enters T v . Let x be the entry vertex. 
Then Q[u, x]T v [x, v] is a u-to-v shortest path in R that does not cross Pi, contradicting our 
assumption. □ 

The algorithm: We can now describe our Bellman-Ford algorithm to relax all edges from 
vertices of C\ to vertices of C2. Pseudocode is shown in Figure HI 

Assume that Rp t and Rp r and distances between pairs of boundary vertices in these graphs 
have been precomputed. 

In each iteration j, we relax edges from vertices of Vc 1 to all v G Vc 2 i n Rp e an d in Rp r 
(lines 9 and 10). Lemma [5] implies that this corresponds to relaxing all edges in R from vertices 
of Vc-i to vertices of Vq 2 ■ By the results in Section 14.21 this suffices to show the correctness of 
the algorithm. 

Lemma H] shows that lines 9,10 can each be implemented to run in 0(| Vcri | + I ^€2 1) ti me - 
Thus, each iteration of lines 6-10 takes ©((IVcJ + \ Vc 2 D^G^Ci I + \ Vc 2 1)) time, as requested. 

It remains to show that Rp r and Rp e and distances between boundary vertices in these 
graphs can be precomputed in 0(\ Vp\ log |Vr|) time 

Shortest path tree T in R with source r\ can be found in 0(|Vr| log | Vp\) time with Dijkstra 
using the feasible price function eft obtained from the recursively computed distances in R. 
Given T, we can find its rightmost path in 0(|Vr|) time by starting at the root r\. When 
entering a vertex v using the edge uv, leave that vertex on the edge that comes after vu in 
counterclockwise order. Computing Rp r given P r also takes 0(|Vr|) time. We can next apply 
Klein's algorithm [5] to compute distances between all pairs of boundary vertices in Rp r in 
0(|Vr| log \Vr\) time (here, we use the non-negative edge lengths in R defined by the reduced 
cost function induced by <fi). We similarly compute Pn and pairwise distances between boundary 
vertices in Rp r 

We can now state our result. 

Theorem 1. Given a planar directed graph G with real edge lengths and no negative cycles and 
given a source vertex s, we can find SSSP distances in G with source s in 0(nlog 2 nj log log n) 
time and linear space. 

Proof. We gave the bound on running time above. To bound the space, first note that finding 
an n/p-division of G using the algorithm of [2] requires 0(n) space. Klein's algorithm [5] and 
Dijkstra also has linear space requirement. The recursively computed distances take up a total 
of 0(p^) = 0{n) space. In the intra-region boundary distances stage, the total memory spent 

on storing distances is 0{p{sj n/p) 2 ) = 0{n). 



In the single-source inter-region boundary distances stage, we need to bound the space for 
our Bellman-Ford variant. The size of each table is 0(b) = 0(n). Since we only need to keep 
tables from the current and previous iteration in memory, Bellman-Ford uses 0(n) space. It is 
easy to see that the last two stages use 0(n) space. Hence the entire algorithm has linear space 
requirement. □ 

5 Concluding Remarks 

We gave a linear space algorithm for single-source shortest path distances in a planar directed 
graph with arbitrary real edge lengths and no negative cycles. Running time is O (n log 2 n / log log n) , 
an improvement of a previous bound by a factor of log log??,. As corollaries, bipartite pla- 
nar perfect matching, feasible flow, and feasible circulation in planar graphs can be solved in 
0(n log 2 nj log log n) time. 

It would be interesting to consider other types of graphs. Results from [1] seem to imply 
that our algorithm generalizes to bounded genus graphs. 

Finding the true complexity of the problem remains open since there is still a gap between 
our upper bound and the linear lower bound. Is 0(n log n) time achievable? 
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Appendix 



The proofs left out or only sketched in the main paper were very similar to those in [6]. For 
completeness, we give them here. 

Proof of Lemma [2] 

We need to show that after iteration j of the algorithm in Figure[TJ ej [v] is the length of a shortest 
path in G from r to v that can be decomposed into at most j subpaths P = P1P2P3 ■ ■ ■ Pj, where 
the endpoints of each subpath Pi are boundary vertices and Pi is a shortest path in a region of 
TZ. 

The proof is by induction on j > (and very similar to the proof of Lemma 4.2 in [6]). 
When j = 0, ej[r] = and ej[v] = 00 for all v € B\ {r} after line 3 and the base case holds. 

Suppose j > and that the lemma holds for j ' — 1. Consider a shortest path P in G from 
r to a v £ B that can be decomposed into subpaths P1P2P3 ■ ■ ■ Pj, where the endpoints of each 
subpath Pi are boundary vertices and Pi is a shortest path in a region of TZ. We need to show 
that after iteration j, ej[v] is the length of P. 

Subpath P' = -P1-P2 • • • jPj—i is a shortest path in G from r to a w £ P> which can be 
decomposed into at most j — 1 subpaths as above. Furthermore, there is a region R £ TZ such 
that v and to are boundary vertices of R and Pj is a shortest path in it! from w to v. 

At some point in iteration j, we reach line 7 with C being the cycle defining the boundary 
of R and v,w £ Vc 1 - By the induction hypothesis, ej_i[V| is the length of P'. Since ej[v] is set 
to a value of at most ej-i[w] + d R (w,v), ej[v] is at most the length of P. 

Let us show the other inequality. For any w £ P>, ej—i[w] is clearly the length of some path 
in G from r to w that can be decomposed into at most j — 1 subpaths, where each subpath is 
a shortest path in a region between two boundary vertices of that region. Hence, when ej [v] is 
updated in line 7, its value is the length of some path in G from r to v that can be decomposed 
into at most j such subpaths. This shows that ej[v] is at least the length of P, completing the 
proof. 

Proof of Lemma [4] 

We need to show that relaxing all edges from Vbi to Vc 2 in Rp can be done in 0(\ Vc x \ + | Vc 2 1) 
time in any iteration of Bellman- Ford. 

We only sketched a proof in the main paper. Let paths P\ and P2 and |Pi| x | P^ | matrix 
A be defined as in the proof sketch. We need to show that the column-minima of A can be 
found in 0(|VcJ + l^cU) time. As shown in the main paper, this amounts to showing that for 
1 < k < k' < \P\\ and 1 < I < V < f(i), we have A ki + A k i V > A kV + A k >i. 

Since P\P P^P is a cycle and since Rp is planar, any pair of paths in Rp from u' k to v[ 
and from u' k , to v[, must intersect in some w £ Vp p , see Figure El Let b k = ej-i[u k ] and 
b' k = ej-i[u k r] (recall that we identified u k with u' k and u k / with u' k ,). Then 

A ki + A wv = (b k + d Rp {u' k ,w) + d Rp {w,v[)) + (b k > + d Rp (u' k ,,w) + d Rp {w,v[,)) 
= (h + d Rp (u' k ,w) + d Rp (w, v[,)) + (b k > + d Rp (u' k ,,w) + d Rp (w, v[)) 
> (h + d Rp (u' h ,v' v )) + (b k > + d Rp {u' kll v[)) 
= A kV + A k i h 



as requested. 




Figure 5: The situation in the proof of Lemma [H Any pair of paths in Rp from u' k to v', and 
from u' k , to v[, must intersect in some w E Vr p . 



